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Abstract 

We review the status of heavy flavour production in QCD. Comparison of 
experimental and theoretical results for top and bottom production are given. 
Selected topics in charm production are also discussed. 



Invited talk given at the 
XV International conference "Physics in Collisions 
Cracow, Poland, June 8-10, 1995. 

CERN-TH/95-256 
IFUM 519/FT 
September 1995 



1 On leave of absence from INFN, Sezione di Milano, Milan, Italy. 
2 On leave of absence from INFN, Sezione di Genova, Genoa, Italy. 



-1- 



Although the first next-to-leading-order calculations of heavy flavour production 
were performed more than five years ago |l 0, progress in this field is constantly 
being made. The work of ref. [[| has confirmed the results of ref. |T|. A calculation of 
the next-to-leading cross section for the photoproduction of heavy quarks has been 
given in ref. ||], and has been confirmed by ref. ||. The computation of the radiative 
corrections to the electroproduction of heavy quarks, via an off-shell photon, was 
presented in ref. 0. A method that accounts for the correlation of heavy quarks at 
next-to-leading order was developed in ref. J7| for heavy quark hadroproduction. An 
application of this calculation to fixed-target production of heavy quarks is given in 
ref. ||. The method was extended to the photoproduction of heavy quarks in ref. 0. 

In the following I will describe some recent experimental and theoretical progress in 
top, bottom and charm production. I briefly remind the reader of the QCD mechanism 
for heavy quark production. The hadroproduction process is depicted in fig. [l| The 




Figure 1: The heavy quark production process in perturbative QCD. 

high energy colliding hadrons can be viewed as a broad band beam of partons (quarks 
and gluons), which collide and fuse to produce the heavy quark pair. Thus the cross 
section is given by the formula 

a QQ = J d %idx 2 f t Hl (xi,/i)/f 2 (x 2 ,/i) (tijfapi, x 2 p2,m Q ,fi,as) (1) 

where a is the short-distance cross section, which is calculable order by order in 
perturbation theory. In fig. [I] only one lowest order diagram is represented. Diagrams 
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like the one of fig. p], after squaring the amplitude to obtain a cross section, give 
rise to a contribution of order a|. Thus, the leading process is of order af , and the 
next-to-leading one is a|. The coupling constant a$ is evaluated at a scale of the 
order of the mass of the heavy flavour. Thus, for charm the coupling constant a$ is 
roughly 0.3 to 0.5, for bottom is around 0.2, and for top is 0.1. We expect therefore 
that theoretical predictions should be very reliable for top, less for bottom, and even 
less for charm. In fact, next-to-leading corrections are around 30-40% for top, 100% 
for bottom, and even larger and less controllable for charm. It is clear therefore that 
for top we expect to be able to predict the cross section rather well. For bottom, and 
expecially for charm, we can expect deviations from theoretical predictions due to 
higher order and non-perturbative effects. In these cases, we also must keep in mind 
whether certain effects have a simple explanation in terms of some non-perturbative 
effects, and whether one can model the cross section with the perturbative result 
supplemented with some model for the non-perturbative effects. 

I will first deal with the issue of comparing the measured top cross section to 
theoretical expectations. I will then discuss the status of bottom production. Charm 
production physics is a rather wide and complex field, mostly because of the fact that 
non-perturbative effects in charm production do play an important role. Thus, there 
are several experimental issues of a certain interest, like the A dependence, and the 
leading particle effect. Here I will focus upon what I think are the most important 
discrepancies between perturbation theory predictions and experimental results. I will 
show that certain distributions, like the transverse momentum of charmed mesons, 
the transverse momentum of the pair, the azimuthal distance of the pair, and the 
invariant mass of the pair, are somewhat in contrast with each other, and that they 
may give an indication of how the comparison of theory and data should be performed. 

1. Top production 

For top hadroproduction, it was found that radiative corrections are generally 
well under control. This allows us to make predictions for top cross sections with a 
relatively small error. The recent measurements of top cross section at the Tevatron 
||10||,|[L1]| have turned out to be in remarkable agreement with the theoretical predic- 
tions. This is illustrated in fig. |2]. In the figure we have reported the result of the 
old calculation of ref. |T2|, using the DFLM structure functions, a next-to-leading 
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calculation using the MRSD- set, and the calculation of ref. |13| which attempts to 



resummation of Sudakov threshold effects. Other calculations are reported in ref [14 



IS |, where in the latter an alternative approach to the inclusion of threshold effects is 



attempted. When estimating the cross section, one should not forget that the uncer- 
tainty on A5, which defines the strenght of the QCD coupling constant, is larger than 
what is usually assumed in structure function fits. In fact, all LEP determinations 



|]16 | tend to give a larger value. Therefore, for comparison, we also show the computed 
cross section in case one uses a value of A 5 more compatible with LEP determination. 
The corresponding increase of the cross section is presumably an overestimate, since 
a larger value of A 5 implies stronger evolution of the parton densities, which soften 
more rapidly as the momentum scale increases. The figure shows that only minor 
changes in the theoretical predictions have taken place in the last few years. We can 
expect that in spite of the fact that more subleading contributions may be added in 
the future to the top strong cross section, they will not change substantially the pure 
next-to-leading result. 




Figure 2: Cross sections for top production: various calculations versus the CDF 
and DO measurements. 
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2. Bottom production 

In bottom production radiative corrections are large, and various estimates of 
corrections of even higher order (as given for example by the renormalization and fac- 
torization scale dependence) lead to theoretical uncertainties of the order of a factor 
of 2. These uncertainties, when combined with other physical uncertainties, such as 
the error in the knowledge of Aqcd and of the structure functions, result in a rather 
poor theoretical prediction. The CDF [17], DO QTSJ and UAl [B| experiments have 



all measurements of the spectrum of B mesons. The earliest CDF measurements 
reported a cross section which was much higher than QCD prediction, and seemed 
to be in contrast with CDF measurements. This problem was due to the poor the- 
oretical understanding of the direct J/^ production, and after the introduction of 
microvertexing techniques, the cross sections have come down to smaller values. The 
remaining discrepancy which is often quoted in CDF publications is due to the fact 
that modern sets of structure functions tend to favour small values for A5, and there- 
fore smaller cross sections. If we instead allow for larger values of A 5 , as favoured 
by LEP experiments, all we can conclude is that all data (including UAl data) are 
consistent with the theoretical prediction, although on the high side of the theoretical 
band. What once seemed to be a discrepancy between CDF and UAl data, is now 
gone, both data sets bearing the same relation with respect to the theoretical curves. 
This is illustrated in fig. ^ The main band (solid) is obtained using the MRSA 
parton densities, except that for the upper line the value of A 5 is taken to be 300 
MeV, not quite consistent with DIS data, but closer to what is indicated by LEP 
measurements. The same curves are also shown without fragmentation effects, in 
comparison with the full band (dotted lines) obtained interpolating table 6 of ref . 0] . 
As one can see not much has changed since then in the theoretical prediction, and the 
only thing one can say is that the CDF measurement is closer to the upper limit of the 
theoretical band. For comparison, I also show data from the DO collaboration in fig. [|. 
The old UAl data is shown in fig. [5[ When comparing theoretical and experimental 
curves, one should remember that certain experimental results are deconvoluted from 
the heavy quark fragmentation effects, and are given as a bare quark cross section, 
while others are heavy flavour meson's cross sections. The former should be compared 
to the theoretical prediction obtained without the inclusion of fragmentation effects, 
while for the latter one should include fragmentation. The UAl experiment presented 
both a B hadron cross section and a h quark cross section, by deconvoluting the effect 
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of fragmentation. In the figure we show both the predicitons of ref. [||, and a more 
modern one, performed using the MRSA structure functions, and the same variation 
of parameters we used at Tevatron energy. Again, we see that the theoretical band 
has not changed much. The UAl quark data is on the high side of the theoretical 
band, as the CDF and DO data. This fact has been quantified in the study of ref. [^(| . 
where it is found that the ratio of the data over theoretical upper band is 0.97 for 
DO, 1.3 for CDF, and 0.84 for UAl. 

We can conclude by saying that data on inclusive b production at hadron colliders 
is consistent with QCD expectation, and lies on the upper extreme of the prediction 
band. In other words, in order to fit the cross section, one is forced to use values of 
the parameters (like the scales, the mass, and the value of A) that favour higher cross 
section. 




Figure 3: Transverse momentum spectrum of B mesons at Tevatron energies, as 
measured by the CDF collaboration. The main prediction band (solid lines) is also 
shown without the inclusion of fragmentation effects. The dotted lines are the old 
prediction of ref. , using the DFLM structure functions sets. 
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As first pointed out in ref. jL|, the perturbative computation of the b cross section 
at hadron colliders reaches a difficult kinematical regime when going from the SppS 
energies to the Tevatron. This is due to the appearance of large logarithms of S/m 2 
(where S is the cm. energy squared, and m is the mass of the heavy quark being 
produced) in the perturbative expansion. In other words, at very large S the expan- 
sion parameter becomes as \og(S/m 2 ), which (for S large enough) becomes of order 
1, thereby spoiling the convergence of the perturbative expansion. Quite substantial 
theoretical work was performed on the small- a; problem [21]. In particular, in ref. [22| 
it was found that the resummation of small-x effects increases the total rate by 30%. 
These effects are very likely to operate in the small p t region. It is unlikely, however, 
that they would strongly improve the agreement between theory and data. 

More detailed analysis of bottom distributions have been performed at the Teva- 
tron. In refs. [|23|] the azimuthal correlation has been measured and compared with 
the theoretical prediction of ref. [0. While in qualitative agreement with QCD predic- 
tion, important differences are observed. We can hope that correlation measurements 




Figure 4: Cross section for b quark production as a function of a transverse momen- 
tum cut, as measured by the DO collaboration. 
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of this kind may shed a light also on the problems observed in the single inclusive px 
distribution. 



3. Charm production 



Figure |^ gives an instructive picture of the uncertainties in charm and bottom 
cross sections at fixed-target experiments. Observe the considerable improvement 
that takes place when going from charm to bottom. Observe also the strong mass 
dependence of the charm result. Needless to say, similar uncertainties plague the 
pN cross sections, as shown in fig. [7]. The experimental measurement of bottom 
production cross section in proton-nucleon collisions at fixed target is a new result 
||24|| . New-generation fixed-target experiments have accumulated very large statistics 
of charm events. For a (possibly incomplete) list of recently published results, see 



refs. [EH] and p6[. For recent reviews of the current experimental situation and future 
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Figure 5: Transverse momentum spectrum of B mesons at 630 GeV center of mass 
energy, measured by the UA1 collaboration. 
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perspectives, see refs. |p7fl , PSfl . For an older review see p9| . Results on bottom 



production at fixed target are given in refs. [3C ] — 131] . As one can see, experimental 
results on total cross sections for charm and bottom production at fixed target are 
in reasonable agreement with theoretical expectations. We remind the reader that 
many puzzling ISR results in pp collisions at 62 GeV remain difficult to explain (see 
the review f2l|), in particular the large production rates reported in ref. pTf . 

Photoproduction results are also in fairly good agreement with theory, as can be 
seen from fig. || The large band in the figure is obtained by varying all parameters, 
including a variation of the charm quark mass from 1.2 to 1.8 GeV. The dashed 
and dotted bands are instead obtained by varying all other parameters, and keeping 
the mass of the charm quark fixed to 1.5 GeV. Very recently a new data point has 
been added by HERA ||32|| . At the Bruxelles conference a new HI result has been 
announced [[33[ . 

Differential distributions for charm are at present in a more complex situation. It 




Figure 6: Cross sections for b and c production in ttN collisions versus experimental 
results. 
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turns out that in photoproduction one finds a remarkably good agreement between 
theoretical expectations and experiments. One typically looks at the transverse mo- 
mentum distribution of a single quark, the transverse momentum of the pair, the 
invariant mass and the azimuthal correlation of the pair. All this quantities are in 
good agreement with theoretical expectations, provided one includes in the calcula- 
tion the effect of a fragmentation function, parametrized in the same way as in charm 
production in e + e~ collisions. This is not the case in hadroproduction. Typically one 
sees that the inclusive p t distribution of a single meson is well described by pure QCD, 
without the inclusion of a fragmentation function. One also sees a similar behaviour 
for the Xf distributions. The azimuthal correlation of charmed pairs requires the 
inclusion of some non-perturbative effects, that could be described as a primordial 
transverse momentum of the incoming quarks. This transverse momentum turns out 
to be of reasonable size, that is to say, below 1 GeV 2 . Another observable which is 
strongly sensitive to a primordial transverse momentum of the incoming quarks is 
the transverse momentum distribution of the charmed pair. In this case, only a very 




Figure 7: Cross sections for b and c production in pN collisions versus experimental 
results. 
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large primordial transverse momentum could reproduce the measured cross section. 
In the following I will illustrate this problems. It is fair to say that a satisfactory 
answer to these problems is not known yet, but also that present data gives some 
hints for a direction in which to deepen theoretical work. I begin by showing in fig. || 
the azimuthal distance between the charm and anticharm, taken from ref. [|34|| . It 
is quite obvious that a strong enhancement in the back-to-back region is observed, 
as expected from a hard production mechanism. Perturbative calculations display a 
similar behaviour, although they are generally much more peaked, as shown in fig. |10[ 
The solid histogram in the figure is obtain by a purely perturbative calculation. The 
dashed and dotted histogram are obtained by assuming that incoming partons do 
have some non-vanishing transverse momentum, which is transmitted to the final 
quark-antiquark pair. This transverse momentum can be viewed as a simple-minded 
model for non-perturbative effects, due to the fact that hadrons have a finite trans- 
verse size. The figure was obtained for a beam energy of 230 GeV, but for the sake 
of the following discussion it can be considered as energy independent. We see that 




Figure 8: Cross sections for c production in 'jN and 7p collisions versus experimental 
results. 
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all data seems to favour a (k^) not larger than 1 GeV 2 , the WA92 data favouring 
even smaller values. A value around 1 GeV 2 is roughly what one would expect, based 
upon experience with Drell-Yan pair production, which requires intrinsic transverse 
momenta of the order of 600 MeV 2 . Drell-Yan pairs are mainly produced via quark- 
antiquark annihilation, while heavy flavours are mostly produced by gluons. Since 
the gluons couple more strongly than quarks by a factor of 9/4, one expects that the 
associated transverse momentum should be somewhat larger. Observe that the az- 
imuthal correlation is independent of fragmentation effets, since fragmentation does 
not change the direction of the outgoing particle. In fig. pi data from the E687 pho- 
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Figure 9: Azimuthal correlation in the production of charmed pairs. 
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toproduction experiment is also reported f35 |. Comparisons of the E687 data with 
theoretical expectation are given in ref. ||, and a good agreement is found, whether 
or not one adds an intrinsic momentum kick to the incoming parton. In fact, in the 
case of photoproduction, a transverse momentum kick has a less dramatic effect, since 
only one parton is carrying the kick, and also because the perturbative distribution 
is already broad, and it washes out the effect of the kick. 

Let us now turn to the inclusive transverse momentum distributions. Here we 
find a problem in the comparison of photoproduction and hadropoduction data, as 
illustrated in figs. |TT] and |T^. We see that photoproduction data is in good agreement 
with QCD predictions, including the effects of fragmentation. Transverse momentum 
kicks do not alter the distribution in an appreciable way. In the case of hadropro- 
duction, one finds instead that the perturbative prediction is in good agreement with 
data, and that the inclusion of fragmentation effects spoil the agreement. In the 
figure, also a transverse momentum kick is included, which improves the agreement, 
but we still see that even a 1 GeV 2 primordial transverse momentum is not enough 
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Figure 10: Theoretical predictions for the azimuthal correlation in the production 
of charmed pairs. 
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to overcome the effect of fragmentation. 

A similar situation arises in the distributions of the transverse momentum of the 
charmed pair. As an example I report in table [I] a computation of various average 
quantities, in pion-nucleon collisions at 230 GeV, in comparison with the recent study 
of ref. |P6 |. These results deserve a comment. We see that, aside from the rapidity 
difference, both the transverse momentum of the pair and the invariant mass of the 
pair would be in better agreement with the data if the effect of fragmentation was 
much less dramatic. In fact, if we don't include fragmentation effects at all, we see 
that all transverse distributions are in good agreement with theoretical prediction, 
provided one is willing to accept the possibility of a non-perturbative intrinsic trans- 
verse momentum effect of the incoming parton of the order of 1 GeV 2 , which is a 
very reasonable value. It is often claimed that the fact that longitudinal distributions 
(like the xf distributions) are harder than theoretical predictions may require the 
inclusion of higher twist effects p7 |. On the other hand, the problem is there only 
if one insists in using the fragmentation functions for longitudinal distributions, a 
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Figure 11: Single inclusive p\ distribution measured by the E769 experiment, com- 
pared to theoretical expectations. 
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procedure that is not fully justified. Many models are capable of justifying easily the 
observed longitudinal distributions. On the contrary, the fact that transverse mo- 
mentum distributions are harder than expected is much more worrying, since in this 
case the factorization theorem should apply. One possible way out of this problem 
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Figure 12: Single inclusive p\ distribution measured by the E687 experiment, com- 
pared to theoretical expectations. 





bare 


<#> = 
no fragm. 


= 1 GeV 2 
with fragm. 


(k 2 T ) = 
no fragm. 


= 2 GeV 2 
with fragm. 


data 




0.277 


2.24 


1.116 


4.2 


2.04 


1.98 ±0.11 ±0.09 


(Mqq) 


4.20 


4.20 


3.62 


4.20 


3.63 


4.45 ±0.03 ±0.13 


(^Vqq) 


0.931 


0.886 


0.669 


0.852 


0.650 


0.54 ±0.02 ±0.24 



Table 1: Various average quantities computed with and without fragmentation 
effects, and transverse momentum kick of 1 and 2 GeV 2 , compared with the results 
of the data analysis of ref. [36 . 
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may be related to the fact that a large fraction of bottom cross section in hadropro- 
duction comes from gluon fragmentation into heavy quarks, while in the case of e + e~ 
annihilation, and also in the case of photoproduction, this fraction is much smaller. 
This possibility requires further study, but it could also explain why collider data for 
bottom production is on the high side of the prediction band. We have in fact seen 
that also in that case, if the fragmentation function is not included, the agreement 
between theory and data becomes much better. 
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